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Abstract 

We discuss a many-particle quantum system, which is obtained by adding some nonher- 
mitian but PT (i.e. combined parity and time reversal) invariant interaction to the A^-i 
rational Calogero model without confining potential. This model gives rise to scattering 
states with continuous real spectrum. The scattering phase shift is determined through 
the exchange statistics parameter. We find that, unlike the case of usual Calogero model, 
the exclusion and exchange statistics parameter differ from each other in the presence of 
PT invariant interaction. 
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1 Introduction 



The subject of exactly solvable many particle quantum mechanical systems with long- 
range interactions is closely connected to diverse subjects like fractional statistics, random 
matrix theory, level statistics for disordered systems, Yangian algebra etc. and generated 
a lot of interest in past years. The A N _i Calogero model (related to A N _i Lie algebra) 
is the simplest example of such a dynamical model, describing N particles on a line and 
with Hamiltonian given by [1] 

i N q2 1 ,.2 JV 

where g is the coupling of long-range interaction and u is the coupling of harmonic confin- 
ing interaction. This An_! Calogero model can be solved exactly to obtain the complete 
set of discrete energy eigenvalues and corresponding bound state eigenfunctions. The 
complete set of energy eigenvalues can be written as 

E ni ,n 2 ,-n N = ^[1 + (N- 1H + U £ Uj , (1.2) 

A 3=1 

where njs are non-negative integer valued quantum numbers with n 3 - < rij+i and v is a 
real positive parameter related to the coupling constant of long range interaction as 

g = v 2 - v. (1.3) 

This spectrum in (1.2) is same as that for a N number of free bosonic oscillators apart 
from a constant shift for all energy levels. This spectrum can be expressed exactly in the 
form of the energy eigenvalues of free bosonic oscillators : E nij7l2t ... nN = 4p + cuJ^jLi^j, 
where fij = rij+v(j — l) are quasi-excitation numbers. These quasi-excitation numbers are 
no longer integers and they satisfy a modified selection rule given by fij +1 — fij>u. This 
selection rule restricts the difference between the quasi-excitation numbers to be at least 
v apart. As a consequence, the Calogero model (1.1) provides a microscopic realization 
for generalized exclusion statistics (GES) [2] with v representing the corresponding GES 
parameter [3, 4, 5]. 

The Calogero model without confining potential, which is obtained by setting uj = 
in eqn.(l.l), is also studied in Ref. [1]. Unlike the earlier case, the spectrum of this 
model is continuous and only scattering states occur. Due to such scattering, particle 
momenta in a outgoing iV-particle plane wave get rearranged (reversely ordered) in terms 
of the momenta in the incoming plane wave. The corresponding momentum independent 
scattering phase shift is given by 9 SC = 7r ^ Af (^~ 1 ) ; w hich is simply utt times the total 
number of two-body exchanges that is needed for rearranging iV particles in the reverse 
order. Thus it is natural to identify v as the exchange statistics parameter in this case 
[4]. It may be noted that this exchange statistics parameter coincides with the exclusion 
statistics parameter as obtained earlier in the presence of confining potential. 
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In the last few years, theoretical investigations on different nonhermitian Hamiltonians 
have received a major boost because many such systems, whenever they are invariant 
under combined parity and time reversal (PT) transformation, lead to either real (when 
PT symmetry is unbroken) or pairs of complex conjugate energy eigenvalues (when PT 
symmetry is spontaneously broken) [6, 7, 8, 9, 10]. Such property of energy eigenvalues 
in nonhermitian PT invariant systems can be related to the pseudo-hermiticity [9] or 
anti-unitary symmetry [10] of the corresponding Hamiltonians. As concrete examples of 
PT symmetric quantum mechanics, the Hamiltonians of only one particle in one space 
dimension have mostly been considered in the literature so far. However, nonhermitian 
but PT invariant extension of some exactly solvable many particle quantum mechanical 
system in one space dimension have also been considered recently [11, 12, 13, 14]. The 
PT transformation for such iV-particle system can be written as 

% — > —i, Xj — > —Xj, Pj — > Pj (1.4) 

where j G [1, 2 • • • N] and Xj (pj = denotes coordinate (momentum) operator of the 

j—th particle. In particular, an extension of A^^i Calogero model with confining term 
is proposed by adding a momentum dependent long-range interaction (5J2 

m _i_ _g-) 

to the Hamiltonian (1.1): 



1 N d 2 a 1 1 d uj 2 N 

= -5 g ^ + 5 E + % <^y Wj + y g A . (") 

where 5 is a real parameter [5, 12]. It is shown that this nonhermitian, PT invariant 
model can be solved exactly and within certain range of the related parameters it yields 
a real spectrum 

Nu N 
E mn2 -n N = —[l + (N-l)u + uJ2n J . (1.6) 

/ 3=1 

Here v = v' — 5 and v' is a real positive parameter which is related to the coupling 
constants g and 5 as 

g = v' 2 - v' (1 + 25) . (1.7) 

In analogy with the case of original Calogero model, the energy eigenvalues (1.6) can 
be rewritten exactly in the form of energy spectrum for iV free oscillators: E nin2 ... nN = 
I-tt+uj X/JLi nj, where fij = rij+u(j — l) are quasi excitation numbers satisfying a modified 
selection rule given by n J+1 — fij > v . Consequently, the extended Calogero model (1.5) 
also provides a microscopic realization for GES, where v represents the exclusion statistics 
parameter. 

In this context, one can naturally ask whether the exclusion and exchange statistics 
parameters are same or not in the case of above described nonhermitian, PT invariant 
extension of the Calogero model. This motivates us to find out the exchange statistics 
parameter associated with such extension of Calogero model in the absense of confining 
interaction, whose Hamiltonian is obtained by putting u = in eqn.(1.5). In Sec. 2 of this 



3 



article, we briefly review the construction of scattering eigenstates and phase shift for the 
An~i Calogero model in absence of confining potential. By following the same procedure, 
in Sec. 3 we find out the scattering eigenstates for the extended Calogero model without 
confining potential and also calculate the related phase shift. Using this scattering phase 
shift, we determine the exchange statistics parameter for the extended Calogero model. 
Section 4 is kept for conclusions. 



2 Scattering states of Calogero model without con- 
fining interaction 

In this section we briefly discuss the scattering state solutions and corresponding phase 
shift of A/v-i Calogero model in the absence of confining interaction. Such a Calogero 
model is described by the Hamiltonian given in (1.1) with u = as 

^ 1 ^ d 2 g ^ 1 ^ ^ 

2 j =1 dxj 2 ._^ k [Xj — Xk) 2 

Following [1], the eigenvalue problem for the above Hamiltonian can be solved to obtain 
scattering states within a sector of configuration space corresponding to a definite ordering 
of particles like x\ > x 2 > ■ ■ ■ > x^. The zero energy ground state wavefunction of this 
model is given by 

w = n(^-^r, ( 2 - 2 ) 

j<k 

where v is a real positive parameter satisfying the relation (1.3). 

Next, we consider the general eigenvalue equation associated with the Hamiltonian 
(2.1): 

W=PV, (2-3) 

where p is real and positive. Solutions of this eigenvalue equation can be written in the 
form 

V> = ip gr T(x 1 ,x 2 , ■ • -x N ) , (2.4) 
where t(xi,X2, ■ ■ -xjv) satisfies the following differential equation, 

_I V — - v V 1 dT = p 2 T (2 5) 

2 j=i dxj [Xj — Xk) dxj 

Further, to separate the 'radial' and 'angular' part of the eigenfunction, one assumes that 
t(xi,x 2 , ■ ■ -xjy) = Pk, q (x)x(r), where the radial variable r is defined as 

r 2 = ^E(x i -x i ) 2 (2.6) 
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and Pk, q (x)s are translationally invariant, symmetric, k-th order homogeneous polynomials 
satisfying the differential equations 

Note that the index q in P^ q {x) can take any integral value ranging from 1 to g(N,k), 
where g(N,k) is the number of independent polynomials which satisfy eqn.(2.7) for a 
given iV and k [1] . One should further note at this point that the translational invariance 
and homogeneity property of the polynomial Pk, q (x) enforce it to satisfy the following 
relations, 

j=l ux 3 j=l UJj 3 

By substituting P^ q {x)x{r) in the place of r(x±, X2, • • • xn) in eqn.(2.5), and using the 
relations (2.7), (2.8), one obtains 



d 2 X (r) (1 + 26) d X (r) 



dr 2 r dr 



= P 2 x(r) , (2.9) 



where b = + k + N{ - N ~ l ^ u , The above equation admits a solution of the form 
x(r) = r~ b J b (pr), where J b {pr) denotes the Bessel function. Hence the scattering state 
eigenfunctions of 7i with eigenvalue p 2 are finally obtained as 

i> = 1] fa - x k Yr- b J b {pr)P Kq {x) . (2.10) 

j<k 

Next we want to discuss the scattering phase shift for the above model. For this 
purpose, one has to construct a more general eigenfunction which in the asymptopic limit 
(i.e., r — > oo limit) can be expressed in terms of an incoming free particle wavefunction 
and an outgoing free particle wavefunction (V>-)> where the incoming wavefunction 
will be of the form 

N 

ip + =exp[iJ2PjXj], (2.H) 
j'=i 

with pj < Pj+i, p 2 = YljLiPj an d X/JLiPj — 0. This can be achieved by taking appropriate 
linear superposition of all degenerate eigenfunctions (with eigenvalue p 2 ) of the form 
(2.10): 

oo g{N,k) 

i^gen = II ( X J ~ x kY E E C kq r~ b J b {pr)P Kq {x) , (2.12) 

j<k k=0 q=l 

where Ck q s are expansion coefficients depending on particle momenta. We assume that 
each of the momenta pi is a product of a radial part p and an angular part ctj, i.e, pi = pcti. 
Now by matching the dimensions of the right hand sides of equations (2.11) and (2.12), it 
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can be shown that C kq = p^~Ck q (cti), where Ck q (ai) depends only on the angular parts 
of the momenta. By using this relation for Ck q and also the asymptotic properties of 
Bessel function at r — > oo limit, 

J b ( pr ) _> _J = / e -<(n+§)f +iPT + e i(„+!)f-iprl 

^271 pr- 
one can write down the asymptotic from of ip gen (2.12) as 

ifjgen ~ ^+ + V- , (2.13) 

where 

oo »(JV,fc) _ 

V>± = (27rr)- V B_3) II - E E C^K^^e^)^ (2.14) 

j<fc fc=0 g=l 

with A = b — k = + N ( N ~ 1 ) U anc l n — ^d^-. Now, to get the expression for the phase 
shift due to scattering, the outgoing wavefunction -0_ has to be expressed in terms of the 
incoming wavefunction For this purpose let us consider a special permutation of the 
particle coordinates T, defined as 

Txi = xjv-i+i , i = 1,2, • • N 

such that the set {— Tx} belongs to the same sector as that of the set {x}, with the 
particle ordering x\ > x<i > ■ ■ ■ > xn- It is also to be noted here that the symmetry and 
the homogeneity of P^ q {x) enables us to write the relation, 

P ktq (-Tx)=e- ilm P ktq (x). 

Hence, by taking advantage of the above facts and finally using eqn.(2.11), it can be shown 
that the outgoing wavefunction ^_ can be written as 



V>_ = e- in( - A+n ^+{x ^ -Tx, p^-p) 



N 

N(N-l) J-L 

= e -^^^e^[iJ2 x jPN+i-j}- (2-15) 

Comparing (2.11) with (2.15), one finds that the momenta of the incoming plane wave gets 
rearranged (reversely ordered) in the scattering process. Moreover, the outgoing plane 
wave acquires a momentum independent phase shift given by -kv n ^ n ~ IS> . Thus v can be 
identified with the exchange statistics parameter associated with this Calogero model. 



3 Scattering phase shift of extended Calogero Model 
without confining interaction 

Here we aim to study the PT symmetric nonhermitian extension of the Calogero model 
without confining interaction. We extend the Hamiltonian (2.1) by adding an extra term 
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b^lj^k (x -- Xk ) af"' which is nonhermitian but symmetric under combined PT transforma- 
tion. Hence the extended Calogero model which we will be considering here is described 
by the Hamiltonian 

1 ^ d 2 1 r ^ 1 d , 

H ex t = -t E JT2 + o E r _ ^ + 6 E 7 _ T o-- ( 2 - 16 ) 

^ j=l OT jr' ^ j^t \ X 3 X k) j-^k \ X 3 Xk > OX 3 

Since H ext (1.5) reduces to the above Hamiltonian at oj = limit, the exchange statistics 
parameter of extended Calogero model can be determined from the phase shift analysis 
of the scattering states associated with the Hamiltonian (2.16). Now we follow the same 
procedure as in the previous section to study the scattering states of this extended model. 

We start by observing that the zero energy ground state wave function of the Hamil- 
tonian (2.16) is quite similar in form with the ground state wave function (2.2) of the 
original Calogero model: 

</v = nfo (2-i7) 

j<k 

where the modified exponent v' is related to the coupling constants g and S through the 
relation (1.7). For the purpose of obtaining nonsingular ground state eigenfunction at 
the limit Xj — > Xj, v' should be a no n- negative exponent. Due to eqn.(1.7), this condition 
restricts the ranges of coupling constants g and S as (i) S > — |, > g > —(6 + |) 2 , and 
(ii) g > with arbitrary value of 5. Next, we consider the general eigenvalue equation 
associated the Hamiltonian (2.16) given by 

n ext ij = p 2 ij, (2.i8) 

where p is a real positive parameter. It is easy to see that the solutions of this eigen- 
value equation can be written in the form = tp gr T'(xi,X2 ■ ■ -xjv), where if) gr represents 
the modified ground state eigenfunction (2.17) and r'(xi, x 2 ■ • ■ xjy) satisfies a differential 
equation like 

1 N B 2 r' 1 Br 1 

_i _ {y> _ 5 ) _ , jf~ = P 2 r' ■ (2.19) 

2 - =1 oXj j_£ k (xj Xk) oxj 

Next we assume that t'(x 1 ,X2 ■ ■ ■ xjv) can be factorised as 

r'(x 1 ,x 2 ---x N ) = P^ q (x)x'(r), (2.20) 

where r is the radial variable defined in (2.6) and (x)s are translationally invariant, 
symmetric, k-th order homogeneous polynomials satisfying the differential equations 

"ff^Jxt l (3 8\ 

g^r +( " -%K^)Wr^) Pt « (x)=0 - (2 ' 21) 

Note that the form of eqn.(2.21) is same as eqn.(2.7) apart from the fact that here v is 
replaced by i/ — S. Hence it is clear that P' k q {x) can be obtained from any given expression 



7 



of Pk t q(x) by simply substituting the parameter v with v' — 5. So the index q in P' kq (x) 
can also take values ranging from 1 to g(N, k). Substituting the factorised form (2.20) of 
t(xi,X2 ■ ■ -xn) in the differential eqn.(2.19) and making use of the properties of P' kq (x) 
we obtain the equation satisfied by the 'radial' part of the wavefunction as 



d 2 X '(r) l + 2Vdx'(r) 



p 2 X\r) (2.22) 



with b' = + k + (i/ - 8) N( - N - 1 \ The solution of eqn.(2.22) can be expressed through 
the Bessel function: x'( r ) = r~ b ' Jy(pr). Hence the scattering state eigenfinctions of H- ex t 
(2.16) with real positive eigenvalue p 2 are obtained as 

1> = II - x k y' r~ v Jv{pr)P^ q {x) . (2.23) 
j<k 

Next, we carry out the calculation of phase shift in the same manner as in the previous 
section. First, we aim to construct a more general eigenfunction of TC ext such that in the 
asymptotic limit it can be expressed in terms of an incoming wave ("0+) of the form 
(2.11) and an outgoing wave (tp-). For this purpose one has to take appropriate linear 
superposition of all degenerate eigenfunctions (with eigenvalue p 2 ) of the form (2.23): 

oo g(N,k) 

= II - *kf E E C' kq r- b ' J bl {pr)P' Kq {x) , (2.24) 

j<k k=0 g=l 

where C' k s are expansion coefficients which are functions of particle monenta. Once again 
by doing dimensional analysis, we obtain 

, (3-N) N(N-1)6 ~ 

C' kq = V 2 2 C' kq (ai) , 

where C' k (ai) depends only on the angular parts of the momenta. By using this explicit 
expression for C' kq and the asymptotic properties of Bessel function at r — > oo, we obtain 
the asymptotic form of ip gen (2.24) as vp gen ~ vp + + ip_, where 

oo g(N,k) 

1> ± = (27TT)" V B '" 5) II - **Y r ~ A ' E E &^)r- k P^(x)e^ + ^^. 

j<k k=0 q=l 

In the above expression A' = b' — k = + (is' — 5) N ( N ~V and n' = + N ( N ~^ 5 , 
Following the same procedure as in the case of Aw-i Calogero model without confining 
term, we find that the outgoing wavefunction (ip~) can be written as 



V>_ = e-^ (A ' +ra 'V+(^ -> -Tx, p -> -p) 



N 

l N(N-l) 



= e- i ™^*-e X p[iT, x jPw-j]- ( 2 - 25 ) 

3=1 
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Comparing eqns.(2.25) and (2.11), we observe that the momenta of the incoming plane 
wave gets rearranged (reversely ordered) in the scattering process. Furthermore, the out- 
going plane wave now acquires a momentum independent phase shift given by ttu' N ^ N 2 ~^ , 
which is tip' times the number of two body exchanges needed for rearranging iV-particles 
in the reverse order. Thus v' can be identified with the exchange statistics parameter 
associated with this extended Calogero model. 



4 Conclusion 



Here we study a nonhermitian but PT invariant extension of the A^-i Calogero model 
without confining interaction. Unlike the case of Ajv_i Calogero model with confining 
term, here we find only scattering states with continuous, real eigenvalues. We also calcu- 
late the scattering phase shift explicitly for this model. The exchange statistics parameter 
for this model is determined through the scattering phase shift of plane waves. Such ex- 
change statistics parameter is given by u' , which is related to the coupling constants 
through eqn.(1.7). However, it is found earlier that the exclusion statistics parameter 
associated with PT invariant extension of Calogero model in the presence of confining 
interaction is given by v = v' — 5. Thus we surprisingly find that, in contrary to the case 
of original Calogero model, the exclusion and exchange statistics parameters differ from 
each other in the presence of PT invariant interaction. In particular, within a range of 
coupling constant given by 5 > 0, > g > —6(1 +5), the value of exclusion statistics pa- 
rameter v becomes negative. As a result, the ground state energy does not have any lower 
bound in the thermodynamic limit (i.e., N — > oo limit). On the other hand, the exchange 
statistics parameter v' is positive and well defined even within the above mentioned range 
of the coupling constants. 
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